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We have systematicaUy studied the hard-core Bose-Hubbard model with correlated hopping on 
a triangular lattice using density-matrix renormalization group method. A rich ground state phase 
diagram is determined. In this phase diagram there is a supersolid phase and a pair superfluid 
phase due to the interplay between the ordinary frustrated boson hopping and an unusual correlated 
hopping. In particular, we flnd that the quantum phase transition between the supersolid phase 
and the pair superfluid phase is continuous. 
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I. INTRODUCTION 

The Bose-Hubbard model, as the first explicit exam- 
ple of quantum phase transition, has been studied ex- 
tensively with many techniques. The simplest version 
of Bose-Hubbard model only involves a Mott insulator 
phase and a superfluid phas&i, and the quantum phase 
transition between these two phases can be well-described 
by a semiclassical Landau-Ginzburg theory. For exam- 
ple, in two dimension, this quantum phase transition is 
either an ordinary 3d XY transition or a z = 2 mean 
field transition depending on the chemical potential. In 
the last few years, it was proposed that various extended 
Bose-Hubbard models can have much richer and more 
exotic behaviors. For example, a Z2 topological liquid 
phase has been discovered in an extended Bose-Hubbard 
model on the Kagome lattice^r-^, and in the same model 
an exotic quantum phase transition between the Z2 hq- 
uid and a conventional superfluid phase was identified -. 
Also, with an extra ring exchange term, it was demon- 
strated both numerically and analytically that the Bose- 
Hubbard model can have an exotic fractionalized Bose 
metal phase^^-S. 

With a hard core constraint, i.e. doubly occupied sites 
are removed from the Hilbert space, the Bose-Hubbard 
model is equivalent to a spin-1/2 model. Due to the rapid 
development of numerical techniques, exotic phases have 
been identified in many quantum spin-1/2 models as well. 
For example, based on the density-matrix renormaliza- 
tion group (DMRG) method, a fully gapped topological 
liquid phase has been discovered in the Kagome lattice 
spin-1/2 Heisenberg modeli^iii, as well as the Ji — J2 
Heisenberg model on the square latticei^. 

In this paper, using the DMRG method, we demon- 
strate that in one simple extended hard core Bose- 



Hubbard model, there are three different interesting phe- 
nomena: first of all, there is a supersolid (SS) phase, 
where there is a coexistence of the off-diagonal long range 
order of boson creation operator, and a boson density 
wave order. Secondly, there is a pair superfluid (PSF) 
phase, where (bi) = while (bibi+a) 7^ 0. This pair su- 
perfluid phase is an analogue of the charge-4e supercon- 
ductor that was discussed lately^iiJ^. Thirdly, we show 
that there is a continuous quantum phase transition be- 
tween the supersolid and the pair superfluid phase. 



II. MODEL HAMILTONIAN 

We consider a hard-core Bose-Hubbard model with a 
correlated hopping on a triangular lattice 

- ^T. 



(uibfbk + h.c. 



(1) 



where bf {bi) is the boson creation (annihilation) oper- 
ator and rii is the boson number operator on site i. In 
this Hamiltonian t is the ordinary nearest-neighbor (NN) 
boson hopping amplitude, and V is the NN repulsive in- 
teraction. i^T is a correlated hopping term, and ijk S A 
are three sites in a small triangle (shown in FigH]) of the 
lattice. We will see that the presence of the K term sig- 
nificantly enriches the phase diagram of this model. In 
the numerical simulations, for simplicity, we will always 
set t = 1 as the unit of energy, and focus on the case with 
V > and K > 0. Notice that here the boson hopping 
term is frustrated. 

We determine the ground-state phase diagram of the 
model Hamiltonian Eq. ([T]) by extensive and highly accu- 
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FIG. 1: (Color online) Illustration of the 3-leg ladder geome- 
try with unit vectors x = (1,0) and y = The boson 
hopping strength is t, and the repulsive interaction is V. The 
correlated hopping term K moves a hard-core boson from one 
site to another one, depending on the number of boson in the 
third site in the same small triangle. 
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FIG. 2: (Color online) Ground state phase diagram for the 
correlated hard-core Bose-Hubbard model (Eq. ([1])) in a tri- 
angular lattice at filling p = 1/6, as determined by accurate 
DMRG calculations on cylinders with Ly up to 9. Chang- 
ing the coupling parameters V and K, four different phases 
are found, including the superfluid (SF), supersolid (SS), pair 
superfluid (PSF), and the phase separation (PS). 



rate DMR G^^'^^ simulations. In particular, we consider 
a system with total number of sites N — ^ Ly, which 
are spanned by multiples L^x and Lyij of the unit vec- 
tors X = (1,0) and y = (i, ^). For the DMRG calcu- 
lation, we consider both a cylinder boundary condition 
(CBC) and a fully periodic boundary condition. Here, 
CBC means open boundary condition along direction, 
while periodic boundary condition along Ly direction. 
This allows us to work on numerous cylinders with much 
larger system size to reduce the finite-size effect for a 
more reliable extrapolation to the thermodynamic limit. 
We keep more than m = 6000 states in each DMRG 
block for most systems, which is found to give excellent 
convergence with tiny truncation errors that can be fully 
neglected. 



III. PHASE DIAGRAM 

The main result of this paper is illustrated in the phase 
diagram of the model (Eq.([T])) at filling p = ^, as shown 
in Figl21 obtained by extensive DMRG studies on nu- 
merous cylinders with Ly = 3 — 9. The nature of the 
ground state of the model Hamiltonian Eq.(IT]) at half 
filling has already been studied previously without con- 
sidering the correlated hopping term K, where a super- 
solid phase was found to be stable over a wide range 
of interaction strength^-"—. In the SS phase there is a 
long range order of both the boson creation operator and 
the boson density wave. Further study^° shows that this 
SS phase will also survive at low boson density, such as 
p — ^, when the repulsive interaction is strong enough. 
For weak interaction, this SS phase will give way to the 
simple atomic superfluid phase. 

In our current paper, we demonstrate that the unfrus- 
trated correlated hopping term will compete with the re- 
pulsive interaction and lead to an interesting phase dia- 
gram. In particular, with intermediate strength of K, the 
SS phase is driven into a uniform pair superfluid phase 
where {bi) — 0, while {bibi+a) 7^ 0. 

To analyze the ground state properties of the system, 
we calculate both the density structure factor 

and the momentum distribution function 

M,(k) = l^e^k.(r.-r,)^^+^^.^^ (3) 

where p is the filling factor of the system. We also calcu- 
late the pair superfluid structure factor 

to characterize the pair superfluid phase. Here Af = 
bibi+a is the nearest-neighbor pair annihilation operator 
along a direction, with a = x, y or y — x. 

In both the SF phase and supersolid phase, the ob- 
tained S{k) and Mii{k) show Bragg peaks at the corners 
of the hexagonal Brillouin zone, e.g., at ki = (±47r/3, 0). 
In particular, at small K, both the peak of the density 
structure factor S{k) and momentum distribution func- 
tion M(k) are very sharp (not shown). As shown in Fig|31 
at large V = 6.0, with the increase of K , S{ki) decreases 
continuously, while Mi,{ki) increases at first and then de- 
creases with larger K. Eventually both S(ki) and Mb(ki) 
becomes very weak beyond certain critical value of K. On 
the other hand, M^(k) shows peaks at zero momentum, 
i.e., fco — (0,0), and will increase with K, Mp{k) in- 
creases monotonically before phase separation. One can 
obtain the corresponding order parameters in the ther- 
modynamic limit based on the finite-size scaling of the 
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FIG. 3: (color online) (a) The density structure factor order 
parameter rria, (b) atomic condensate density Cb at momen- 
tum fci = (47r/3,0), and (c) pair condensate density Cp at 
momentum fco = (0,0), as functions oi K ?A, V = 6.0 and 
filling p = 1/6 , respectively, with system size A'^ = 12 x 6, 
12 X 9, 18 X 9, and the corresponding extrapolations in the 
thermodynamic limit. 



FIG. 4: (color online) Examples of finite-size scaling of (a) the 
density structure factor order parameter ms and (b) atomic 
condensate density Ch at momentum fci = (47r/3, 0), as well as 
(c) the pair condensate density Cp at momentum fco = (0,0), 
for different K a,t V — 6.0 and filling p — 1/6, with system 
size up to 18 X 9. 



peak values 5(fci), Mf,(fci) and Mp{ko). Specifically, the 
structure factor order parameter and the boson conden- 
sate density can be determined by rris = S{ki)/N and 
Cb — Mi,{ki)/N, while the total pair condensate density 
is given by Cp = E„ = Ea M^{ko)/N. Examples of 
these order parameters are shown in FigjH] as a function 
of correlated hopping K with — 12 x 6, 12 x 9 and 
18 X 9. 

Nonzero rn^ and Cb in the thermodynamic limit cor- 
respond to the diagonal long-range order (LRO) and off- 
diagonal long-range order (ODLRO), respectively. Ex- 
amples of the finite-size scaling are shown in Fig|4] by 
plotting iris, Cb and Cp as functions of l/N, at = 6.0 
and different K, using quasi-2D lattice geometry with 
system size up to = 18 x 9. The obtained order pa- 
rameters extrapolated to the thermodynamic limit are 
presented in FigE] (dark cyan sphere) and FigH) It is 
worth noting that after extrapolation both rus and Cb 
continuously decrease to zero simultaneously at the criti- 



cal K , and the system enters the pair superfluid phase in 
the thermodynamic limit, indicating a continuous phase 
transition. However, Cp is nonzero in both the SS and 
PSF, and it always increases monotonically with K (See 
Fig. [Sl^c)) before phase separation. Specifically, the 
finite-size scaling for Cp at = 2.0 in the supersolid 
phase gives us a small but finite value Cp « 2.5 x 10"''' 
(around 2% of Cp at K = 4.8), as shown in FigHl^c). 

With the decrease of the repulsive interaction V , both 
the supersolid phase and pair superfiuid phase become 
weaker and move to smaller K. In particular, a super- 
fluid phase appears at small K locating at the left bottom 
corner in the phase diagram (see Fig|2]). Compared with 
the supersolid and pair superfiuid phase at large V , the 
pair condensate density at ko = (0,0) becomes zero, i.e., 
Cp(fco) = 0- Even with the absence of the repulsive in- 
teraction V , the atomic superfluid phase is completely 
suppressed with the increase of K and give way to the 
supersolid phase through a first-order transition around 
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FIG. 5: (color online) (a) The density structure factor order 
parameter rus, (b) atomic condensate density Ct at momen- 
tum fci = (47r/3,0), and (c) the pair condensate density Cp 
at momentum fco = (0,0), as functions of it^ at = 0.0 and 
filling p = 1/6, respectively, with system size A*' = 12 x 6, 
12 X 9 and 18 x 9, as well as the corresponding extrapolations 
in the thermodynamic limit. 



K = 0.3. This can be seen clearly in FiglSJ in which 
TTis and Cp jump from zero to a finite value in the ther- 
modynamic limit, while Cb encounters a sharp drop to a 
smaller but still finite value. It is interesting to note that 
the pair superfluid phase still remains robust in a finite 
parameter region K w 1.6—2.0, even in the absence of the 
repulsive interaction V, due to the competition between 
t and K. Finally, the system becomes phase separated 
when K becomes dominant. 



IV. THEORETICAL UNDERSTANDING 

In this section, we develop a simple Landau- Ginzburg 
theory to describe the supersolid, superfluid and paired 
superfluid phases on the triangle lattice, including the 
phase transitions between them. The dispersion of a bo- 
son with frustrated hopping in triangular lattice has two 



minimums at K and —K, where K = (47r/3, 0). There- 
fore we introduce the two U{1) order parameters (/)i.2 for 
boson condensate at ztK 



b^x) ^ M^y"""" + M^)e- 



iK-x 



Under time reversal and lattice symmetries, 
form as follows: 



(5) 

hi 2 trans- 



e 

uii) 



h 4>2, 02 01 

h 02, 02 4>1 

h 02, 02 01 

bi 016*'^, 02 026' 



i2e*""/3, z = 1,2,3 



(6) 



Here Ti are translations by the three unit vector a.; of the 
triangle lattice: ai = (1,0), a2 = (—1/2, a/3/2), = 
(— 1/2, — -\/3/2). Cq is 7r/3 rotation. is reflection 

about X axis. Q is time-reversal transformation. C/(l) is 
the global phase transformation. 

The supersolid phase is characterized by two nonzero 
condensate order parameters 0i and 02, which have mo- 
mentum K and —K respectively. The coexistence of 
these two superfluid order necessarily induces a nonzero 
density wave order p ^ 0i02 at momentum 2K = —K, 
where we have used the fact 3K = up to a recipro- 
cal lattice vector in triangular lattice. The presence of 
both superfluid and density-wave order taken together 
signals a supersolid phase, as we found numerically. On 
the other hand, the paired superfluid phase is character- 
ized by a nonzero pair condensate order parameter A at 
zero momentum, whereas both 0i and 02 are disordered. 

Since the supersolid phase has a lower symmetry than 
the paired superfluid phase, the phase transition can be 
understood as the development of long-range single bo- 
son superfluid in the paired superfluid (disordered) phase 
within the framework of Landau-Ginzburg theory in 2-f 1 
dimension. The effective Landau-Ginzburg Lagrangian, 
whose form is dictated by the symmetry property (6), is 
given by 
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As r decreases, 0i and 02 become nonzero, which signals 
the onset of superfluid order as well as the associated 
density wave order 0J02- Note that the pair order pa- 
rameter A is nonzero across the transition. Without loss 
of generality, A is chosen to be real and positive. Due to 
the trilinear coupling term between A, 0i and 02, C is 
minimized by 02 = 0i = 0, where is complex. The last 
term in S is symmetry allowed for the triangular lattice, 
and locks the relative phase between 0i and 02 , and pins 
the phase of p = 0i02 = 0^ to three distinct values cor- 
responding to three degenerate density wave patterns in 
the supersolid phase. In terms of 0, £ is given by 



|a^0|2-fri|0|2-Hu'|0r+^'('/'' + 



(7) 
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Except for the last term, L is the standard complex scalar 
field theory in 2+1 dimension, and this transition belongs 
to the 3d XY universality class with order parameter 0. 
By power-counting, the sixth-order phase-locking term is 
strongly irrelevant at the critical point, thus this transi- 
tion is continuous. 

Interestingly, the above XY transition differs from con- 
ventional paired to single boson superfluid transition, 
which lies in the Ising universality class. The distinction 
arises from the fact that the superfluid phase studied here 
has two coexisting (rather than one) condensate order 
parameters, of which the relative phase is a well-defined 
physical quantity associated with the density wave order. 
Our theory can be directly tested by further numerical 
studies on the critical exponent for the physical order pa- 
rameter at the quantum critical point. For example, the 
critical exponent v take the value of the ordinary 3d XY 
transition: v ^ 0.67. Also, the density wave order pa- 
rameter p '-^ 0^ is a bilinear of in Eq. [71 thus our theory 
predicts that the order parameter p has an anomalous di- 
mension r\ ~ 1.49. These predictions can be verified by 
further numerical studies on model Eq. [1] 

We can also interpret the supersolid to pair superfluid 
transition in terms of the topological defects inside the 
supersolid phase. For convenience, let us rewrite the 
fields 01 and introduced in Eq. [5]as follows: 

(/-l = ^, 02 = 0* ^- (8) 

and 0* are complex fields that carry lattice momentum 
K and —K respectively, while -0 carries the U(l) sym- 
metry of the original boson operator hi. Thus the boson 
density wave order parameter is p ^ 0i02 ^ s-nd the 
pair superfiuid order parameter is bibi^a ^ 0i02 ^ 0^- 

In Eq.m because 0i and 02 are both physical degrees of 
freedom, and ip are defined up to a Z2 gauge ambiguity, 
i.e. the physics is unchanged under transformation — > 
— 0, "0 — >■ —ip. The supersolid phase corresponds to the 
case where both and 4' are condensed, while in the 
pair superfluid phase only ^ is condensed. Inside the 
supersolid phase, the smallest superfluid vortex, has only 
TT— vorticity, i.e. it is a bound state between a tt— vortex 
of tp and a tt— vortex of 0. In other words, both tp and 
will change sign after encircling this vortex, while the 
physical degree of freedom is unchanged. Notice that the 
TT— vortex of is a full vortex of boson density wave order 
p, which is equivalent to a dislocation of the density wave 
pattern. 

Starting with the supersolid phase, if we want to drive 
a transition into the pair superfluid phase, we need to 



"melt" the boson density wave order by condensing its 
defects. However, since the pair superfluid phase also 
has a half quantum tt— vortex, this transition cannot be 
driven by condensing the smallest tt— vortex discussed in 
the previous paragraph. Instead, it must be driven by the 
condensation of the 27r— vortex of 0, which only melts the 
boson density wave, but leaves the superfluid stiffness un- 
affected. Since the physical boson density wave order pa- 
rameter p ~ 0^ is a bilinear of 0, this transition is driven 
by condensing a "double" dislocation of p. Thus this 
transition is analogous to the transition between pair- 
density-wave and charge-4e superconductor discussed in 
Ref^. 



V. SUMMARY AND EXTENSION 

We have dstudied a hard-core Bose-Hubbard model 
with an unusual correlated hopping on a triangular lat- 
tice using density-matrix renormalization group method. 
In the phase diagram, we discovered a supersolid phase 
and a pair superfluid phase, in addition to the stan- 
dard superfluid phase. The supersolid and pair superfluid 
phases were discussed separately before in different spin 
models^"— 1^. However, to our knowledge it is the flrst 
time that both these phases are realized in one model. We 
also theoretically establish that the phase transition be- 
tween the supersolid phase and the pair superfluid phase 
is continuous. 

If the phase coherence and superfluid stiffness of the 
pair superfluid phase are destroyed, then the system most 
likely enters a fully gapped Z2 liquid phase with the same 
topological order as the toric code model^i. Presumably 
this new transition can be obtained by turning on some 
extra terms in the Hamiltonian. We will leave this to 
future study. 
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